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ABSTRACT 

Models of large scale (magnetohydrodynamic) dynamos (LSD) which couple large 
scale field growth to total magnetic helicity evolution best predict the saturation of 
LSDs seen in simulations. For the simplest so called "a 2 " LSDs in periodic boxes, 
the electromotive force driving LSD growth depends on the difference between the 
time-integrated kinetic and current helicity associated with fluctuations. When the 
system is helically kinetically forced (KF), the growth of the large scale helical field is 
accompanied by growth of small scale magnetic (and current) helicity which ultimately 
quench the LSD. Here, using both simulations and theory, we study the complementary 
magnetically forced(MF) case in which the system is forced with an electric field 
that supplies magnetic helicity. For this MF case, the kinetic helicity becomes the 
back-reactor that saturates the LSD. Simulations of both MF and KF cases can be 
approximately modeled with the same equations of magnetic helicity evolution, but 
with complementary initial conditions. A key difference between KF and MF cases is 
that the helical large scale field in the MF case grows with the same sign of injected 
magnetic helicity, whereas the large and small scale magnetic helicities grow with 
opposite sign for the KF case. The MF case can arise even when the thermal pressure 
is approximately smaller than the magnetic pressure, and requires only that helical 
small scale magnetic fluctuations dominate helical velocity fluctuations in LSD driving. 
We suggest that LSDs in accretion discs and Babcock models of the solar dynamo are 
actually MF LSDs. 



1 INTRODUCTION 



Understanding the magnetohydrodynamics of large scale 
magnetic field generation in turbulent astrophysical rota- 
tors is the enterprise of large scale dynamo (LSD) theory. 
The LSD problem is usually posed in kinetically forced (KF) 
circumstances, whereby an initially weak field is subject to 
relatively strong hydrodynamic forcing. LSDs then describe 
the growth or sustenance of fields on spatial or time scales 
large compared to the largest scale of the underlying tur- 
bulent forcing. Small scale dynamos (SSD), in contrast, de- 
scribe field growth on scales at or below the turbulent forcing 
scale. LSDs and SSDs are contemporaneous and interacting. 
Understanding that interaction, and how KF LSDs saturate 
have been the subject of much research. 

The LSDs of class i c 20th c entur y textbooks on the 



subjec t (lMoffattMl978l) . iParkerl (|l979l L iKrause fc Raedlerl 



l|l980h . lZeTdovichl fl983r0 do not predict LSD saturation. A 
related fact is that these original LSDs do not conserve mag- 
netic helicity, the inclusion of which has proven effective at 
improving the prediction of LSD saturation. Progress in un- 
derstanding principles of KF LSD saturation has emerged 
from simple numerical experiments to compare with theo- 
retical predictions of mean field models. The simplest ex- 
periments are those of the so called a 2 dynamo in which 
isotropically driven kinetic helicity is injected into a peri- 



odic box at intermedi ate wave number k = 5 and the k = 1 
large scale field grows (|Brandenburd (|200lf q The saturation 
of the LSD in such simulations is best modeled by theories 
that solve for the coupled evolution of large scale helical 
magnetic field and the small scale helical magnetic field. In 
such models, the electromotive force (EMF, £ ) of the mean 
field theory emerges to be the difference between kinetic he- 
licity and current helicity, each multiplied by a correspond- 
ing corr elation time. This difference was first proposed spec- 
trally in lFrisch et alj (|l975l ) and derived and solved dynam- 
ically as a two-scale mean field theory in Black man fc Field! 
(2002). The results of the dynamical mean fiel d theory agree 
with t he saturation seen in the simulation of iBrandenburd 
|200l| y Recently the three scale version of the theory was 
also c ompared to numerical simulations (|Park fc Blackmanl 
(2012). Further extensions and approaches for understand- 
ing LSDs via tracking magnetic helicity evolution and flow 
for more realistic open systems that include boundary flux 
terms is a fruitful ongoing enterprize. 

The basic insight gained from the basic periodic box KF 
a 2 LSD studies is that as the kinetic helicity is forced and 
drives large scale helical magnetic fields, the small scale heli- 
cal field grows because magnetic helicity is largely conserved. 
The build up of the small scale current helicity associated 
with the small scale field then off-sets the driving kinetic he- 
licity and quenches the LSD. This physical mechanism and 
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the equations that describe it also motivate consideration of 
the complementary case of driving the initial system with 
current helicity rather t han kinetic helicity. Thi s was in- 
vestigated analytically in Blackman fc Field! \2004 ) where it 
was found that indeed driving with current helicity produces 
a magnetically forced(MF) analogue to the a 2 dynamo. In 
this paper we present simulations of this MF analogue to 
the a 2 dynamo and compare the results with a two-scale 
mean field theory. We will force the system in the induc- 
tion equation with an electric field that drives small scale 
magnetic helicity. Because the MF system is driven mag- 
netically, there is no "kinematic" regime in the sense that 
the magnetic field is strong from the outset. However the 
initial velocity fluctuations are weak. Thus the analogue of 
the kinematic regime for the MF a 2 case is a "magnematic" 
regime where the small scale velocity fluctuations are small. 
As we will see, both simulations and theory show that it 
is indeed the buildup of the kinetic helicity that ends the 
magnematic regime and quenches the MF analogue of the 
a 2 dynamo. 

Simulations of an inverse cas c ade o f magnetic helicity 
was also found in lAlexakis et alJ (|2006l ) but with a focus 
on non-local transfer functions in the interpretation of the 
inverse cascade, and without a dynamical theoretical model 
in terms of a mean field theory. Here we will present new 
simulations, and compare the simulations and theory. We 
also address the astrophysical relevance of MF LSDs. 

In section 2 we more specifically discuss the basic prob- 
lem to be solved and the computational methods. In section 

3 we present the basic results of the simulations. In section 

4 we develop the theory in more detail and discuss the com- 
parison between theory and simulations. In section 4 we dis- 
cuss why MF LSDs are in fact of basic conceptual relevance 
to LSDs of accretion discs, Babcock type stellar dynamos, 
and magnetic relaxation of astrophysical coronae and lab- 
oratory fusion plasma configurations. We also discuss some 
basic open questions. We conclude in section 5. 



2 PROBLEM TO BE STUDIED AND 
METHODS 

We used the hi gh-ord er finite difference Pencil 
Code (Brandenburg] (|200ll )l along with the message 
passing interface(MPI) for parallelization. The equations 
solved by the code are 

Dp „ 

■m = -' v - u ' 

Du 2 J x B 

-7^7 = ~c s VlogpH 

Dt p 



dA 
~dt 



+ ^(V 2 u+ iVV-u) 
P 3 

ux B-?7J + f. 



(1) 

(2) 
(3) 



where p is the density, u is the velocity, B is the magnetic 
field, J is the current density, D/Dt = d/dt + u ■ V, is 
the advective derivative, L rf is magnetic diffusivity, 'v' is 
viscosity, and c„ is the sound speed, and f is forcing function. 
The default units of the pencil code are 'cgs'. 

We employ a periodic box of dimensionless spatial vol- 
ume (2-7r) 3 and a mesh size of 216 3 . No large scale velocity 



forcing is imposed. Instead, our forcing function / is placed 
in the magnetic induction equation and is a fully helical, 
Gaussian random force-free function(V x f oc f) given by 
f(x,t) = ATfk(t) exp [i k(i) ■ x + icj)(t)], where TV is a nor- 
malization factor and k(t) is the forcing wave number with 
|k(t)| ~ 5 ( 4.5 < |k(t)| < 5.5) and |f k |=0.03(|f k | of 0.07 
was used in l|Park fc Blackman! {{201$ ))). This forcing func- 
tion can be thought of as imposing a small scale external 
electric field e ext such that f = e ex t- It is a source term 
for magnetic helicity for the small scale field proportional to 
{e ex t ■ b). 

From the solutions of the equations, we compute a range 
of quantities including the box averaged spectra of kinetic 
energy(_Bfe), magnetic energy(SM), kinetic helicity((v ■ V x 
v)) and magnetic helicity((A ■ B)). From these quantities 
we further construct mean quantities for comparisons with a 
two-scale theoretical model as discussed in subsequent sec- 
tions. We note that although our system is magnetically 
driven, the "plasma beta" (i.e. /3 P = ratio of magnetic en- 
ergy density to thermal energy density) satisfies P P > 1 most 
of the time in our simulations. 

We will consider Prandtl numbers of unity in this pa- 
per, for two different values of the resistivity. In some of 
the analysis of our equations we need to integrate over all 
wavenumbers and it will be noted that in both cases we 
take kmax to be 107. For the two values of r\ that we choose, 
namely T) = 0.0025 and f] = 0.006 we find that k dis ~ 107 
for rj — 0.006 and k max < kdis for r\ = 0.0025 but that very 
little power in any of the relevant quantities appears above 
k = 107 for either case. Thus we typically used the same 
upper bound in k for both cases. 

For the data analysis, we used Interactive Data Lan- 
guage(IDL). 



3 SIMULATION RESULTS 

3.1 Basic Results: Growth of Large Scale Field 
and Saturation 



Fig [T(a)p(b)1 show a comparison between our simulation re- 
sults and two-scale theory (discussed further in section 3) 
for the dimensionless large and small scale current helici- 
ties (kfhi and k%h,2 respectively), kinetic energy (e), kinetic 
helicity (h v ), and EMF (= £ = (v X b) ~ Q). The simu- 
lation data for the large scale field correspond to quantities 
at ki = 1 and the simulation data for the small scale field 
correspond to the sum of quantities for 2 < k < k m ax- The 
analytic two scale model (discussed later in more detail) 
takes ki = 1 for the large scale and the forcing scale ki — 5 
for the small scale. 

Roughly, the large scale current helicity growth in in 
Fig. 1 re sembles that of the curre nt helicity studied in the 
KF case (|Park fc Blackman! (|2012[ )). Unlike the KF case, the 
the large scale magnetic (and current) helicity has the same 
sign as that on the small scale for the present MF case. This 
is because our forcing in the induction equation injects mag- 
netic helicity of a fixed sign into the system and magnetic 
helicity is conserved until resistive terms kick in. The MF a 2 
LSD thus acts as a dynamical relaxation process: the injec- 
tion of small magnetic helicity drives the system away from 
equilibrium and the LSD competes with this injection by 
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Figure 1. (a) and (b): Time evolution of dimensionless current helicity of large scale field from simulations (blue curve) and small scale 
field (purple curve) compared with theoretical predictions (dotted lines) form our 2-scale model for the large scale k\h\ and small scale 
k^h2 current helicities. The first two panels are for the two different magnetic diffusivities shown. Current helicity is normalized in units 
of b^/k2(Hi = hib^/k2) and time in units of k2b r (t = rk2b r ). (c) Large scale k^h\ from panels (a) and (b) shown on the same plot. 
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Figure 2. (a) and (b) Dimensionless kinetic energy e = {v 2 )/{b 2 ), kinetic helicity h v = (v ■ V X v)/(k2{b 2 .)) and electromotive force 
(EMF) : Q = £ii/(6^), for two different magnetic diffusivities. 



converting much of the small scale injected magnetic helic- 
ity into large scale magnetic helicity of the same sign. The 
magnetic energy of a fixed amount of magnetic helicity is 
less on large scales than on small scales and much of the 
energy lost in the relaxation goes into kinetic energy. Note 



from Fig 2(a) and 2(b) that a sizeable fraction of the kinetic 



energy gained is helical. 

As we will see in the next subsection, most of kinetic he- 
licity is on the forcing and smaller scale. The growing kinetic 
helicity ultimate quenches the LSD in the non-linear regime. 
This mechanism of quenching is supported by a comparison 



between k\h,2 curves in Fig 1(a) and 1(b) with h v curves in 
Fig |2(a)| and |2(b)j In dimensionless form, the LSD of section 
4 is driven by the difference between h,2 and h v . When these 
two are roughly equal (the difference being only equal to 
smaller diffusion and resistive terms) significant quenching 
is expected. Dividing the \t2h2 curves by k\ = 25 shows in- 
deed that the value /12 ~ h v in the late time regimes of the 
plots. In this context, note that Fig lc shows that the higher 
Reynolds number (equal to magnetic Reynolds in this pa- 
per) case takes longer for the large scale field to saturate. 
This is analogous to the KF case where the higher magnetic 
Reynolds number case also takes longer to saturate. However 



for the MF case, it is the viscosity not the resistivity that 
damps the back reaction because the back reactor is kinetic 
helicity. In contrast, the KF case saturates because the small 
scale current helicity builds up, which is damped by resistiv- 
ity. The analytic dependenc e on magnetic Prandtl nu mber 
for the MF was discussed in iBlackman fc Field! (|2004l ). We 
do not pursue that further here. 



3.2 



Time Evolution of Energy and Helicity 
Spectra 



For the MF case, Figs 3(a)|4(c) show the time evolution of 



the magnetic helicity, magnetic energy and current helicity 
spectra as the simulation for the two different values of diffu- 
sivities used. Note the strong early peak at the forcing scale 
and the ultimate transfers to the k = 1 scale. 

Figs 5(a)|6(by show the corresponding spectral evolu- 
tion of the kinetic helicity and kinetic energy. Note that the 
kinetic energy is less peaked that than the magnetic energy 
at the forcing scale for this MF case and that there is some 
growth of kinetic energy at the large k — 1 scale. As we 
discuss later, this is due to the Lorentz force at k — 1 before 
the large scale field becomes force free. 



4 Kiwan Park and E. G. Blackman 



Fig 7(a) and 7(b) compare the kinetic and magnetic en- 
ergy spectra on the same plots and Figure 8(a) shows the 
time evolution of the kinetic and magnetic energy spectra 
for the MF case of the present pap e r com pared with the 
KF simulations of IPark fc Blackmanl (|2012h . Note that the 
KF case shows a stronger peak at the forcing scale for the 
kinetic energy than the MF case, and the MF case shows a 
stronger peak at the forcing scale in magnetic energy than 
the KF case. 

To study the spectral evolution of various quantities 
and their migration between scales, it is useful to define 
average wave numbers with the different bas es defined below 
l|Davidsonl (120041 ). IPark fc Blackmanl <|2012h ) as seen below: 



energy is dominated by that on the large scale. Fig 13(c) and 
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where ki will be taken either as ki = 1 or ki = 2 depending 
on whether we compute the averages for all wave numbers 
or extract the large scale k = 1 in order to compute small 
scale averages. 

Analogously, for (fe) v .v, (k) v . u , and (fc) w . w 
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The time evolution of the average wave numbers in 
the kinetic bases is shown in Fig |9(a)1 and Fig |9(b)1 for 
the MF case, which can be compared with the KF case 
(|Park fc Blackmanl (|2012T l) of Fig |10(a)| and Fig |10(b)| The 
time evolution of the average wave numbers in the magnetic 



bases for the MF case is shown in Fig 11(a) an d Fig| 11(b) 
which can be compared with KF case in Fig |12(a)| and 
Fig 12(b) I. Each of the figures shows two kinds of averages 
the "total average" for which the integration range extends 
from ki = 1 to k max = 107 or the "small scale" for which 
the integration range extends from ki — 2 to k max = 107. 
To distinguish these, we use the subscript "tot" to indicate 
averages defined by the former and "s" by the latter. For 
example, we write fcj-B, tot or fcj.B.s, or ktot and k s when not 
specifying the basis. 

Using k = 1 for the large scale and the range k = [2, 107] 
for the small scale, Fig 13(a) and 13(b) highlight the increase 
of small scale magnetic energy and current helicity at early 
times that results from the MF case followed by the longer 
term growth of the large scale field. Eventually the magnetic 



Fig 13(d) show the analogous curve for the kinetic energy 
and kinetic helicity. Although the kinetic energy and kinetic 
helicity also grow in response to the MF, the kinetic energy 
is dominated by small scale contributions. 



Fig 5(b) and 6(b) show that most of the kinetic energy 
remains on small scales during the simulations but there is 
some growth of large scale kinetic energy. There is no in- 
dependent invers e cascade of velocity in 3-D, only in 2-D 
(|David son (2004)) so the growth of any large scale veloc- 
ity in our case is expected to be the direct result of the 



Lorentz force on those scales. Indeed according to Fig 9(a 



9(b) we can see that the average wavenumbers k v 2 tot and 
k v 2 3 are unequal, highlighting that there is some energy in 
the large scale k = 1 velocity. This inequality disappears as 
time passes and coincidence occurs approximately when the 



magnetic field becomes fully helical(Fig 14(a) 14(b) 15(a) 



15(b) I. Since the magnetic field can only transfer energy to 
velocity fields when the Lorentz force(J x B, Eq. (|18[0 is fi- 
nite on this scale, no energy can be transferred to the kinetic 
eddy once the field becomes fully helical. The evolution of 
the helicity fractions of both the kinetic and magnetic he- 
licity fractions as a function of wave number and time are 



shown in Figs 14(a) and 14(b) The time evolution of the he- 
licity fractions for k = 1 are shown in Fig 15(a) and 15(b) I 



The evolution of the mean wavenumbers seen in 



Fig 11(a) and 11(b) highlights the evolution of the magnetic 
energy and helicity spectra in mf case. At very early times 

fcj-B, tot = k B 2 tot — fcAB.tot = fcj.B.a = k B 2 s = fcAB,s = 5, 

i.e., ktot^ks in all bases. This means the magnetic energy 
resides almost exclusively at the forcing scale(fc/) in this 
early time regime. From Eq.@, the average k at this state 
is simply fcA • B/A ■ B|fc = fc,, which is '5'. Since magnetic 
field and current helicity are helical(J = fcB = fc 2 A), the 
basis independent profile in the early time is explained. 

By t ~ 50, the energy begins to cascade both for- 
ward (to small scales) and inversely (to large scales). For 
~ 50 < t <~ 100, we still have fctot~fc s in all bases (since 
all 6 curves of of e.g. Fig 11(a) I continue to overlap in this 
time regime) but (fc)A-B, (fc)B B, and (fc)j.B split. The fact 
that the quantities all decrease together in this regime shows 
that although there magnetic energy is still overwhelmingly 
contained in the small scales, the quantities are evolving to- 
ward larger scales but they have not crossed over to the large 
scale yet. The two regimes just discussed can be referred to 
as "pre-LSD I" and "pre-LSD II" regimes because there is 
little growth of magnetic energy at fc = 1. 

Beyond t ~ 100 Fig |ll(a)| and |ll(b)| show the influ- 
ence of the LSD. Note that the evolution of k s and ktot for 
each quantity now deviate from each other. That is, calcu- 
lating the fc averages is dramatically affected by changing 
the lower integration bound from fc = 1 to fc = 2. By this 
time, b, v, and EMF have grown and driven growth of the 
large scale k — 1 field. Also at this time, the sign of dk s /dt 
changes from negative to positive. For example, (fc)A B,s of 



T] = v = 0.0025 (Fig 11(b) |, transitions from negative to 
positive growth over the regime ~ 100 < t <~ 200. The 
reason for this change is two fold: (i) The LSD takes in- 
versely transferred magnetic helicity from the small scales 
to fc = 1, thus dropping magnetic helicity and magnetic en- 
ergy out of the bins that contribute to fc^fEg fTTj) . and (ii) 
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Figure 3. Time evolution for the spectra of (a) magnetic helicity (b) magnetic energy (c) current helicity for the r\ = 0.006 run. 
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Figure 4. Time evolution for the spectra of (a) magnetic helicity (b) magnetic energy (c) current helicity for the 77 = 0.0025 run. 



The small helicity is continuously driven at the forcing scale 
k = 5. 

The buildup of kinetic helicity eventually slows the 
LSD. Already around t ~ 200 the slopes of the curves in 
Fig |ll(aJ| and |ll(b)| flatten, seemingly consistent with the 
time at which (v • V x v) approaches saturation and approx- 
imately equals (j • b) . That marks the end of the kinematic 
regime, as discussed in the previous subsection. Eventually, 
the system reaches a steady-state. 



3.2.1 Relation between ktot and k s 

We can mathematically relate ktot and k 3 for each basis. For 
example, for the case of /ca b, tot and &a.b,s (we let E = 
^fcmax f Qr conven i ence j n this subsection), we have 

_5>ab + lAB E feab 



given by 



k A -: 



Eab+AB ' " E a b 

=> k A -B, s ^ a ■ b + A ■ B = k A -T3,tot ( ^ a ■ b + A ■ B 
=>• (Aa-b,s - kA-B,tot) ^2 a b = (k A -B,tot - 1)^4 ■ B 
A • B fc A -B,s — k A . Bi tot 



E< 



k A -I 



1 



(10) 



When fcA B,s=fcA B,tot, the numerator on the right hand 
side(RHS) is 0. This means most of the magnetic helicity 
is in the small scale. In contrast, when fcA B.tot ~ 1, most of 
the magnetic helicity is in the large scale. Generally the time 
rate of change of the large scale over small scale helicities is 



d_ 

dt 



A B 



(fcAB,s — kAB,tot)(kAB,tot — 1) 



(&AB, tot 
(feA-B,s 



I) 2 

kA-S.tot) ftA-B.tot 



(&A-] 



(11) 



In the early time regimes I and II, /ca b.s = ^a b, tot 
and kAB,tot = A;ab,s, so 'd (A • B/E a ' b)/dt' is zero. 
Once /sa b.s / kA B, tot and kA-B,tot < feA B.s), this quantity 
becomes positive, i.e., the inverse cascade sets in. 

The equations of this subsection can also be applied to 
kj.B,tot and fcj.B.s, and can be used to quantify the shape 
of the curves in e.g. Fig | 11(a)] and [Tl(b)| For example, when 
kA B,tot < and /ca-b.s = 0, the inverse cascade begins to 
accelerate. 



4 THEORETICAL TWO SCALE MODEL 
4.1 Basic Equations 

The equations in this pap er were modified from the theory in 
iBlackman fc Fiel d (2004) with respect to the external force 
term. In that paper, the forcing was imposed such that the 
small scale magnetic helicity was kept constant from t = 0. 
In the present case, to better match what was done in our 
simulations, we impose the force in the induction equation. 
The modified magnetic induction equation becomes 



^ = Vx(VxB - // J ~ F 



(12) 
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Figure 5. Time evolution for the spectra of (a) kinetic helicity and (b) kinetic energy for the r) = 0.006 
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Figure 6. Time evolution for the spectra of (a) kinetic helicity and (b) kinetic energy for the 



: 0.006 run. 




where F is the forcing function. We divide the system into 
large scale (wave number fci, or fez,) and small scale(wave 
number fe, or k s ). 

We use an overbar to indicate the mean and lower case 
letters to indicate fluctuations. We assume that the forcing 



applies only to the small scale so that F = f (i.e. F = 0), and 
is helical such that | V x f | — | fca f" | We also assume the mean 
velocity is zero such that V = V + v = v. We can then write 
the remaining MHD quantities as E = E + e, B = B + b, 
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Figu re 8. Time evoluti on of the E^ n and E mag spectra for the MF case of r\ = 0.0025 plotted together with that of the KF case 
iPark fc Blackman|[2012T) . 



J = J + j. The electric field (E = -V x B + 77 J - F) in large 
and small scales are 



E = -f + ryJ, 
and 

e = -v x B + £ + j] j — f , 



(13) 



(14) 



where £ = v x b. 

The magnetic helicities in the large and small scale are 



dt 
and 

dHj 
dt 



= -2(E ■ B) = 2{£ ■ B) - 2??(J ■ B) 



(15) 



= -2(e ■ b) ~ -2{£ ■ B) + 2r/(j ■ b) + 2<b ■ f). (16) 



0£ 



The equation for £\\{=£ ■ B/|B|) can be obtained using 
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which i s derived analogously to that in l|Blackman fc Fieldl 
l|l999h . IBlackman fc Fieldl (|2002T )') from the small scale mo- 
mentum and magnetic induction equations given by 



9v — 

— = v x (V x v) - v x (V x v) - VP e/ / + j x B 



-J x B + vV v, 



(18) 



and 

.au _ 

-pt = B-Vv-v-VB + b-Vv-Vx£: + AV 2 b + k 2 f . 
at 

(19) 

The equation for the EMF £ becomes, 



b • V x b — v-Vxv 



B 

W\ 



-c ki v 2 £ 11 + ki v x f ■ . 

IBI 



3 IBI 



(20) 



1/2 

Here the term involving c v 2 £\\ is the result of the "min- 
imal tau" closure, used e.g iBlackman fc Fieldl |2002t ) that 
incorporates triple correlations T' as well as microscopic 
magnetic diffusivity '77', kinetic viscosit y V, t riple and the 
last term (Eq.(6) in IBlackman fc Fieldl (|2002t )'). The quan- 
tity c is a constant that will later be determined empirically. 

The mean kinetic energy per unit mass (v 2 /2= (w 2 )/2) can 
be obtained from the momentum equation, 

d ^ V "> 2(v ■ j x B) + 2(v ■ J x b) + 2v(-v ■ V 2 v) 



dt 



2((v x j) ■ B) - 2((v x b) ■ J) - 2i/((Vv) 2 )(21) 



Using vector iden t ities and defining the helicity ratios 
<|Park fc Blackmanl (|2012D ). 



mi — 



fh 



|B,I 2 



v,-Vxv, 



J. B, 

ki I B^ 1 2 



(i = 1,2), 



hv 2 

we find the equation for {v 2 } to be 
d(v 2 ) 



dt 



2{f m2 k 2 - f m iki)(S ■ B) - 2vk 2 2 {v 2 ). 



(22) 



(23) 



(24) 



We nondimensionlize the equation with the scalings be- 



low: 



H 1 = hiH 2s , H 2 = h 2 H 2sy r = tk 



3/2 



'H 2 



Q = £ n /k 2 H 2s ,e = U/k 2 H 2s ,R M = M2 > 

R v = Hl{ 2 /vky\H 2s = b 2 rms /k 2 . (25) 



The use of k 2 and root mean squared small scale magnetic 
field b rrns {= br) in these normalizations is natural since the 
external MF was applied to k — k 2 — 5, and the RMS 
small scale field is steadily injected. The ordinary differential 
equations to be simultaneously solved are the dimensionless 
versions of Eq. (fT5j) . ([TBI . ([20)l . and ([2TJ) , which are given 
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Figure 9. Mean kinetic wavenumbcrs in the different kinetic bases for the two different resistivity cases for MF case. Simulations used 
a forcing parameter /o=0.03 
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Figure 10. Mean kinetic wavenumbers in the different kinetic bases for the two different resistivity cases for KF case . Simulations used 
a forcing parameter / =0.07 
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Figure 11. Mean magnetic wavenumbers in the different magnetic bases for the two different resistivity cases for MF case. Simulations 
used a forcing parameter / =0.03 
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Figure 13. Time evolution of dimensionless (a) magnetic energy E m ag, (b) current helicity (A;^ ft;); (c) kinetic energy(e), and 
helicity({v ■ V X v)) in the large and small scale. The large scale (A; = 1) kinetic quantities e and (v ■ V X v) are less than 1% 
the small scale(fc = [2 ~ 107]). 



(d) kinetic 
of those in 



respectively by 

i(£ s 



dhi 
dr 



dh 2 
dr 



dQ 
dr 



and 
de 



_2 

/ml 



2 ( k i\ 1,2 r,uV 2 



2r? kl 

r k'2 



(26) 



A/2 
J ml 



Or 



5 (forcing) 



—h 1/2 f 

,1/2 A 2 ; , 



(27) 



6 (a~r(<b-Vxb>-(v-Vxv>)) 



0,1/2 
J ./ml 



( fe - A2e) (|)"V'-i/i?^) 



3/2 



e h 1 



1/2 




(28) 



10({(vXj)-B)-{(vxb)-J>) 



f l/2 
J ml 



kl \ /fcl 



1/2 ^Q_|^ 2e .(29) 



In our analytic two-scale model we have assumed that the 
kinetic helicity of the fluctuations satisfies (v ■ V x v) = 
fh2k2{v 2 } with fh2 fixed and k 2 = kf = 5. Here the sign of 
the small scale kinetic helicity is positive, which is consistent 
with our simulations expected when H^ 1 is initially positive 
from the imposed forcing. Nondimensionally, we have /12 = 
/h2£. Although the simulations show that fh2 at k — 5 is 
not strictly constant at all times during the simulations, 
the assumption of such is not too far off. In more accurate 
treatments a separate dynamical equation for /12 should be 
included. 

Note also that the term labeled 5 on the right of Eq. (|27|l 
resulted from (b ■ f) and was not solved for dynamically. We 
empirically determined the magnitude of /', the magnitude 
of the forcing needed in the theoretical equations to best 
match the simulations which employed a forcing magnitude 
/o- The magnitude of forcing used in simulation /o = 0.03 
is slightly larger than the values of /' employed in the two- 
scale model (see Table 1). 

The term labeled 9 in Eq. (l2^ namely (v x f) ■ | was 
taken to be zero: it requires a correlation between two nearly 
isotropic functions. This assumption is consistent with the 
simulations as the quantity is measured to be even smaller 
than (v x b), which also requires a correlation between two 
nearly isotropic functions. 
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Figure 14. Comparison of the time evolution of fractional kinetic and fractional magnetic helicity as a function of k for the two different 
Ry = Rm cases. This kinetic energy on the large scale saturates at very low fractional helicities compared to the magnetic case. Note 
also from the previous figure that there is negligible kinetic energy on the large scales compared to magnetic energy; only the fractional 
helicities are plotted in the present figures. 
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(a) (b) 
Figure 15. Time evolution of the fractional helicities fh,kin an d fh,mag on the large scale k = 1. 



4.2 Discussion of Solutions 

We numerically solve the ordinary set of differential equa- 
tions Eq.([26]), (271), (HU), and {29j) for hi, h 2 , Q and e (from 
which we obtain h v ). Table 1 shows the parameters defined 
in the previous section that best match the simulation. The 



dotted lines in Figs 1(a) 1(b) 1(c) 2(a) 2(b) represent the 



solutions with these parameter choices. 

From Eq. (|28|) for Q, we see that the a coefficient (term 
6) increases at first due to h 2 supplied by the forcing. This 
drives the growth of hi of the same sign and the kinetic 
energy e also grows. Some of this kinetic energy is helical 
and so h v grows. This h v grows with the same sign of h 2 so 
when the former is subtracted from the latter in a the overall 
alpha effect is reduced. Thus h v acts as the back reaction in 
a. 

The associated buildup of e also increases the turbu- 
lent diffusion of Q (7) and the dissipation term (8). Note 
that the growth of the turbulent diffusion coefficient as a 
backreaction in the MF case differs from that of the KF 
where the driving itself directly supplies steady turbulent 
diffusion coefficient from the outset. The turbulent diffusion 
has no magnetic component and so must wait for the veloc- 



r\ = v 


V r 


Re M 


/ml 


fm2 


tin 


b r 


0.006 


0.23 


48 


1 


0.615 


0.645 


0.88 


0.0025 


0.28 


141 


1 


0.4 


0.185 


1.5 






c 


/' 


Qo 


hio(10- 14 ) 


h 20 


eo(10" 8 ) 


0.006 


0.45 


0.02 


0.042 


1.23 


0.015 


9.3 


0.0025 


0.8 


0.0195 


0.022 


0.062 


0.0054 


5.7 



Table 1 . Table of parameter values used in the analytic equations 
to best match the simulations. Initial values Qo, hio, and /120 are 
dimensionless values (Eq.j25j). We use ki = 1 and k 2 = 5. 



ity field to build up. This delay in growth of the turbulent 
diffusion is consistent with the difference between the early 
time evolution of the mean magnetic basis wavenumbers in 
ll(a)| and |ll(b)| compared to KF case of 



the MF case for Fig 



Fig 12(a) and 12(b) 



In the latter, the KF immediately leads 
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to a turbulent cascade available to induce scale separation 
in the magnetic quantities. 

As the backreaction from the kinetic energy and kinetic 
helcity grow, dQ/dt converges to zero and Q becomes satu- 
rated at small positive value after the initial bump (Fig 



2 a 



2(b) I The right hand side(RHS) of Eq.([26j (terms 1+2) then 



also converges to zero and h\ saturates. The behavior of /12 
can be explained similarly. 

Complete saturation occurs when all of the terms on 
the right sides of the dimensionless equations Eq. (|26p . (|27[) . 
(|28p . and (|29|) equal zero. However, the largest backreaction 
term is the kinetic helicity term and so setting /12 ~ h v gives 
an approximate assessment of quenching. This approximate 
equality is consistent with the asymptotic values of /12 and 
h v shown in Fig 1(b) and 2(b) given that k 2 = 5. The late 



time growth is limited by viscosity in the sense that the 
faster h v viscously dissipates, the less effective its backreac- 



tion is on /12. This is evident in Fig 1(c) which shows that 
the lower Rv = Rm case approaches saturation earlier than 
the Rv = Rm case. 

It is important to emphasize that it is the viscosity and 
not the resistivity that limits the MF LSD. This contrasts 
the KF case which is resistively limited by how fast the 
small scale magnetic helicity resistively dissipates relative 
to the large scale magnetic helicity. In that respect, keeping 
Rm fixed while reducing Rv increases the magnetic Prandtl 
number a nd the saturation ra t e of t he MF case. This was dis- 
cussed in iBlackman fc Field! (|2004T l. We do not explore the 
Prandtl number dependence in the present paper. In general, 
future work will be needed to better understand the early 
time regime and the associated Rm and Rv dependence. 



4.3 Analytic Simplifications 

Although Eq.(g5J, (gTJl l|25). and (O require a numerical 
solution, some insight can be gained by an analytic approx- 
imation. The hi equation has the form of a Bernoulli equa- 
tion and can be linearized by changing variables. If we use 



111 1 
dui 



l/2x 



the differential equation for hi becomes, 
1 



Tj kl 1 /fciNVa 

or = -^ i+ 7^y q^-p^+ 

J ml 



(30) 



b r &2 ' 



■I ml 

The solution of this standard inhomogeneous differen- 
tial equation is known . Using the same method as in 



IPark fc Blackmanl (|2012r i we find, 



Ul(t) 



- /* p x dt' 



-Pit 



J* Pldt' 



Qi(t')dt' + Const 



(c 



Pit 



,/Qi(*') 



Qi(t) 
Pi 



+ ^Const — 



Pi 

Qi(Q) 
Pi 



dt' + Const 



-Pit 



(31) 



As f -> 00, 'ui(t)(= h\^ 2 y approaches to 'Qi(i)/Pi'(= 
£{i)ky 2 1 f^ibrT/k^ 2 ) which can be calculated by setting 
'dhi/dt' to be zero. This approximate value l Qi(t)/Pi' 



1I1 

matches the numerically solved after t ~ 1500, for 

n = v = 0.006. 

1I1 

The equation for h 2 can not be linearized. Instead, 



setting dhi/dt = 0, and 112 
equation: 



1/2 



we obtain a quadratic 



u 2 {t) 



f 



Vk2b r f m / 2 \ \Vk2brf^2 



f 



Ab r Qh\ /2 fkiy/ 2 
(32) 



V.fmi 



Except the early time regime(i < 130, r\ = v = 0.006), 
this solution fits the numerical result well. These ap- 
proximate equations explain how the saturated magnetic 
helicities in large and small scale depend on '77'. 



5 ASTROPHYSICAL RELEVANCE OF 

MAGNETICALLY FORCED LARGE SCALE 
DYNAMOS 

The MF LSD that we have discussed in this paper main- 
tains a plasma /3 P greater than unity for most of the time. 
Although we consider the simplest such MF LSD in that we 
employ periodic boundaries and no rotation or shear, the 
basic concept of a LSD that is driven by small scale mag- 
netic fluctuation rather than kinetic fluctuation in j3 p > 1 
environments has important conceptual relevance to both 
solar dynamos and dynamos in accretion discs. 

In the solar context, helical LSD models have been sep- 
arated into two classes: (1) helical forcing is primarily ki- 
netic helicity driven by thermal convection (2) flux trans- 
port models wh ich helical fields are the result of magnetic 
instabilities (cf. ICharbonneaul (|2007r 0. We suggest that the 
second class of dynamos can be viewed as an MF LSD, albeit 
with more complexity than the simple version we consider 
in the present paper. Magnet i c inst ability driven dynamos 
in the radiative zone jSpruij \200i )) also seem to fit into 
this category. 

For accretion discs, large scale dynamos are 
now commonly s ee n in s hearing box simul ati ons(e.g. 
iBrandenburg et al.l (tl995h; iDavis et all l|2010l ): [g rcsscl 
I 20ld ); iKapyla, fc Korpil ^Oljy T The magneto-rotational 
instability operating in these simulations produces turbu- 
lence that drives an LSD. But the the magnetic fluctuations 
typically exceed the kinetic fluctuations and it seems 
that the best agreement between mean field theory and 
simulation requires that the driver of the LSD growth 
is not the kinetic h elicity term but the current helicity 
term (|Gressell (|201Cf 0. Thus accretion discs are another 
environment where more general /3 P > 1 MF LSD dynamos 
operate. 

Although the calculations of the present paper are in 
the P P ~> 1 limit, we note that for the /3 P 1 limit, 
MF LSD have long bee n studied in the laboratory context 
fe.g. |ji fc Prage r (2002)). The direct analogy to these low 
f3 MF LSDs may occur in astrophysical coronae (Blackman 

dloo3)). 
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6 CONCLUSION 

We performed numerical simulations of the analogue of an 
a 2 LSD in a periodic box when the system is driven with 
magnetic (or current) helicity rather than kinetic helicity. 
The simulations indeed show that LSD action results from 
the injection of small scale magnetic helicity, analogously to 
the KF LSD action from injection of kinetic helicity. This 
analogy is expected because the growth driver is propor- 
tional to the time integrate difference between current he- 
licity and kinetic helicity. 

We compared the simulation results with a two scale 
theory and found general consistency with respect to the 
basic mechanism of large scale field growth and saturation. 
When the system is magnetically driven at k = 5, the large 
scale k = 1 helical magnetic field grows in the MF with 
the same sign as that on the driving scale. Injecting k — 5 
magnetic helicity drives the system away from its natural 
relaxed state and the LSD evolves the magnetic helicity to 
the large scale where the same amount of magnetic helicity 
has lower energy. Eventually the LSD saturates because the 
growth driver is the difference between kinetic and current 
helicities and the kinetic helicity builds up to quench the 
LSD. This situation complements the KF forced case where 
the LSD drives oppositely signed growth of large and small 
scale magnetic helicities and the current helicity quenches 
the LSD. 

We presented the time evolution of the spectra of mag- 
netic energy, magnetic helicity, current helicity from the sim- 
ulations. Taken together, these spectra exhibit the expected 
inverse cascade of magnetic helicity that is at the heart of 
LSD action, and the absence of an inverse cascade of ki- 
netic helicity. Future work is needed to better understand 
the early time growth regime and study the dependence on 
Reynolds number and Prandtl numbers. 

The MF LSD studied herein was for MHD plasmas with 
ratios of thermal to magnetic pressure overall larger than 
unity. We checked that this was the case at all times during 
the simulations. We discussed that such an MF LSD may 
ultimately be involved in producing the LSD action in MRI 
simulations and may also be away to distinguish Babcock 
type solar dynamo models from KF kinetic helicity driven 
solar dynamo models. 
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